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Abstract 


Puberty is characterized by hormonal, physical and psychological transformation. The human brain undergoes significant changes between 
childhood and adulthood, but little is known about how puberty influences its structural development. Using a longitudinal sample of 711 mag- 
netic resonance imaging scans from 275 individuals aged 7-20 years, we examined how subcortical brain regions change in relation to puberty. 
Our regions of interest included the amygdala, hippocampus and corpus striatum including the nucleus accumbens (NA), caudate, putamen and 
globus pallidus (GP). Pubertal development was significantly related to structural volume in all six regions in both sexes. Pubertal development 
and age had both independent and interactive influences on volume for the amygdala, hippocampus and putamen in both sexes, and the caudate 
in females. There was an interactive puberty-by-age effect on volume for the NA and GP in both sexes, and the caudate in males. These findings 
suggest a significant role for puberty in structural brain development. 


Abbreviations: NA, Nucleus Accumbens; GP, Globus Pallidus 
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Introduction 


The past 15 years have seen a major expansion in research on the structural development of the human adolescent brain, based largely on the 
results of cross-sectional and longitudinal magnetic resonance imaging (MRI) studies (Brain Development Cooperative Group, 2012; Giedd etal, 
1996; Lenroot et al., 2007; Ostby et al., 2009; Raznahan et al,, 2011; Sowell et al., 2002). Studies of brain growth trajectories over adolescence to 
date have predominantly considered growth in relation to chronological age, with few exceptions (Paus et al, 2010; Raznahan et al. 2010). 
However, this approach ignores the cascade of developmental processes during adolescence that are related to puberty, including sexual matu- 
ration, linear (height) growth, body fat re-distribution and maturation of many other physiological systems (Patton and Viner, 2007; Wheeler, 
1991). Since there is a normal variation of 4-5 years in the timing of onset of puberty in healthy humans (Parent et al,,.2003), pubertal develop- 
ment is partially dissociable from chronological age. Examining brain development in relation to pubertal maturation may provide additional in- 
formation regarding the mechanisms associated with adolescent brain development. In this study, we investigated how the developmental trajec- 
tories of subcortical regions that are linked to stereotypical behaviours are associated with pubertal development (Forbes and Dahl, 2010; 
Steinberg, 2008): the amygdala and hippocampus, which play an important role in emotion and mood regulation (Davidson et al, 2002); and the 
corpus striatum including the nucleus accumbens (NA), caudate, putamen and globus pallidus (GP), which are involved in decision-making and 
reward-seeking behaviours (Gottfried,.2011). 


It has been hypothesised that the brain restructuring and development seen in adolescence may be specifically related to the hormonal influ- 
ences that control the onset of and progression through puberty (Giedd et al, 1999; Lenroot et al., 2007; Peper et al. 2011: Sowell et al.,.2002). 
Sex steroids such as testosterone (an androgen) and oestradiol (an oestrogen) have been shown to be capable of inducing both synaptogenesis 
and synaptic pruning in rats and nonhuman primates (Ahmed et al, 2008; Hajszan et al, 2008; Sato et al, 2008) with differential effects of an- 
drogens and oestrogens on different brain areas, which may be related to hormone receptor distribution (Clark et al., 1988; Sholl and Kim, 
1989). These differential effects across brain areas may provide an explanation for the diverging growth trajectories of particular brain struc- 
tures between males and females documented across studies, and the resultant increasing sexual dimorphism in adolescence reported in some 


regions (Brain Development Cooperative Group, 2012; Lenroot et al, 2007; Neufang et al., 2009; Sowell et al. 2002). 


Cross-sectional studies have reported that puberty is associated with aspects of brain development in adolescence. A study focussing on the as- 
sociation between brain volumes and both pubertal stage and testosterone concentration found that males and females in later stages of pu- 
berty, and with higher circulating testosterone concentration, had larger amygdala volumes and smaller hippocampal volumes than their less 
well developed peers (Neufang et al,, 2009). In contrast, a second study investigating puberty and pubertal hormone correlations with grey mat- 
ter volume replicated the positive association between amygdala volume and pubertal stage in boys, but showed decreasing amygdala volume 
with increasing testosterone levels in girls (Bramen et al. 2011). These studies have limited power due to their relatively small sample sizes and 
cross-sectional methods. To date, no longitudinal studies incorporating pubertal measures have been published. Longitudinal analysis allows 
comparison of brain volumes both between-subjects and also within each subject over time, and therefore can provide a measure not just of 
brain volume at a particular time-point, but also of developmental trajectories for each of these subcortical brain regions by following what hap- 
pens to each participant. This is particularly advantageous when looking at brain volumes, which vary substantially between individuals (Brain 
Development Cooperative Group, 2012; Giedd et al, 1999; Tamnes et al. 2013). 


Using a large sample of scans from 275 individuals scanned longitudinally, we examined how subcortical brain regions change in relation to pu- 
berty as measured by Tanner stage (Tanner and Whitehouse, 1976), and compared developmental trajectories in females and males. We used a 
large dataset containing information on pubertal stage and chronological age to examine growth trajectories over adolescence in each of these 
subcortical regions. Based on previous cross-sectional findings, we predicted that the volume of the amygdala and hippocampus would increase 
between 7 and 20 years (Giedd et al., 1999; @stby et al, 2009), whilst the volumes of the corpus striatum - the NA, caudate, putamen and GP - 
would decrease (Brain Development Cooperative Group, 2012; Ostby et al, 2009; Sowell et al,,2002). We hypothesized that volume change for 
all structures of interest would be related to pubertal development as measured by Tanner stage, and that puberty and age would have indepen- 
dent effects on volume in these regions. Therefore, we hypothesized that models incorporating both Tanner stage and chronological age as ex- 
planatory variables for volume change would provide a significant fit for the developmental data from the structures examined. 


Experimental procedures 


Participants 


The sample used for the analysis was taken from the NIMH Child Psychiatry Branch Section on Brain Imaging longitudinal dataset of structural 
MRI scans (Giedd et al., 1996). This large dataset consists of more than 6500 scans from more than 3000 participants of which approximately 
half are typically developing and half are from various diagnostic groups. Inclusion and exclusion criteria for the overall dataset can be found on 
the NIMH website (http://clinicalstudies.info.nih.gov/cgi/detail.cgi?A 1989-M-0006.html). For each participant, the outcomes of interest that 
were measured included ethnicity, socioeconomic status (using Hollingshead scales), IQ (estimated using age-appropriate Wechsler Intelligence 
Scales) and handedness (using Physical and Neurological Examination of Soft Signs inventory (Denckla, 1985)), each of which was collected at 
the time of the first scan, and pubertal status, ascertained at the time of each scan using Tanner stage diagrams (Taylor et al,,2001). Tanner stag- 
ing is the most widely used technique for assessing physical puberty stage. For each participant at each time-point, a single combined Tanner 
stage score was assigned based on the overall stage that the participant felt best described themselves from looking at the separate 
breast/genital and pubic hair scores. 


The subgroup used for the current analysis consisted of 275 unrelated individuals (117 females), and incorporated the scans of all individuals 
who met the following criteria: 


1. Healthy individuals at the time of scanning. Participants were screened for neurological or psychiatric illness using a telephone screening 
interview and completion of a parent-report screening questionnaire (CBCL) at each time-point, and only healthy participants were included 
in the analysis as established by standardized scoring (Achenbach, 1991). All participants included had an IQ of greater than 80. 

2. Two or more MRI scans between the ages of 7.0 and 20.0 years. This age range was estimated using two large US population-based studies 
on pubertal timing (Herman-Giddens et al., 1997; Sun et al., 2002), and incorporating ages from two 2 SD below the mean age of being in 
Tanner stage 2 (6.3-7.8 years for females and 6.9-9.7 years for males depending on the measure (breast/pubic hair/gonadal development) 
used), to 2 SD above the mean age of being in Tanner stage 5 (20.5-21.0 years for females, 19,6-20.1 years for males depending on the 
measure used). This range was refined based on our dataset to 7.0 to 20.0 years, as this incorporated the ages at which there was variation 
in Tanner stage between individuals. 

3. Complete data for age and self-reported Tanner stage for each MRI scan. Individuals who reported regression of puberty during adolescence 
(as indicated by reducing Tanner score) were excluded from the analysis (N = 1; female), since pubertal regression is essentially biologically 
implausible and likely represents error. Pubertal arrest, which may be perceived as regression, is associated with significant systemic illness 
or malnutrition. 

4. Only one individual per family. The original dataset incorporates a high number (> 400 individuals) of monozygous and dizygous twin pairs, 
as well as siblings. The heritability of brain structure and development is still not well-understood, but there are likely to be both genetic and 
environmental effects on structural brain volume and developmental trajectory, as well as pubertal timing, which are not independent 
between family members, and therefore might bias the analysis. Where data were available for more than one family member, only one 
individual was included in the analysis; the family member included was determined by the participant with a higher number of high quality 
scans, more complete demographic data, or, if these were equal, by random selection. 


Of the 275 participants whose data were included in the study, 87.3% were right-handed (7.3% left-handed, 5.4% mixed). The majority (89.5%) 
were Caucasian (5.1% African-American; 2.2% Hispanic; 0.7% Asian; 2.5% other). Details of socioeconomic status, IQ, puberty status and num- 
ber of scans can be seen in Table 1. 


Participants were recruited from the community through local advertisement and reimbursed for their participation in the study. The study re- 
search protocol was approved by the Institutional Review Board of the National Institutes of Health and written informed consent and assent to 
Participate in the study were obtained from parents/adult participants and children respectively. 


Image acquisition 


All MRI scans were T-1 weighted images with contiguous 1.5 mm axial slices and 2.0 mm coronal slices, obtained on the same 1.5-T General 
Electric Signa scanner (Milwaukee, W1) using a 3D spoiled gradient recalled echo sequence with the following parameters: echo time: 5 ms; rep- 
etition time: 24 ms; flip angle: 45°; acquisition matrix: 256 x 192; field of view: 24 cm. The same scanner, hardware and software were used 
throughout the scanning period. All scans were assessed by a clinical neuroradiologist for gross abnormalities. All scans performed as part of 
the NIH project have been rated for motion-related quality by trained technicians. Only scans given a high-quality rating were included in the 
analysis. 


Image processing 


Subcortical volume estimation was performed with the Freesurfer 5.1 image analysis suite using the programme's automated segmentation pro- 
cedure (http://surfernmrmgh.harvard.edu/). This procedure has been described in detail previously (Fisch] et al., 2002), and is summarized 
here, An optimal linear transform is computed that maximises the likelihood of the input image, given an atlas constructed from manually la- 
belled images. A nonlinear transform is then initialized with the linear one, and the image is allowed to further deform to better match the atlas. 
Finally, a Bayesian segmentation procedure is carried out, and the maximum a posteriori estimate of the labelling is computed. The segmenta- 
tion uses three pieces of information to disambiguate labels: (1) the prior probability of a given tissue class occurring at a specific atlas location; 
(2) the likelihood of the image given that tissue class; and (3) the probability of the local spatial configuration of labels given the tissue class. The 
automated segmentations have been found to be statistically indistinguishable from manual labelling (Fisch] et al 2002), and correlations be- 
tween Freesurfer segmentation and manual labelling of hippocampal volume by trained researchers reached 0.85 in a study by Tae et al, (2008). 


Our analysis focussed on the following subcortical brain regions, based on a priori hypotheses regarding changes in volume over adolescence 
and pubertal effects: the amygdala, hippocampus, NA, caudate, putamen and GP. Each of these regions is defined by the automated Freesurfer 
segmentation procedure based on location, likelihood of tissue class and spatial configuration. An example of T1W scans showing both the raw 
data and with Freesurfer's automated segmentation can be seen in the Supplementary information (Fig. S1). 


Statistical analysis 


We conducted our analysis using the average volume across hemispheres, to produce one value for each ROI. Previous studies have demon- 
strated no evidence of developmental difference between hemisphere in these regions (Brain Development Cooperative Group, 2012; Ostby et 
al,,_2009), and our own dataset shows high correlations between hemispheres for all volumes (r = 0.5-0.9, p < 0.001). We analysed raw volumes 
for each region and percentage change for each volume over time. Mixed effects modelling was used (R version 3.1-102; nlme package (Pinheiro. 
etal,,2011)) to analyse the data, thereby allowing an estimation of the fixed effects of measured variables on volume change, whilst incorporat- 
ing the longitudinal nature of the data by including within-person variation as nested random effects. Age was centred on 7 years, which repre- 
sented the minimum age included in the sample. Tanner stage was treated as a continuous variable for this analysis, allowing the model to incor- 
porate changing Tanner stage and changing brain volume for each individual. 


For each structure, volume was first modelled against Tanner stage (Tanner-only model), and linear, quadratic or cubic developmental trajecto- 
ries were modelled. Males and females were modelled separately, to allow for different trajectories of growth through adolescence (Lenroot et 
al. 2007). Tanner stage was treated as a continuous variable to allow the model to account for movement of individuals between stages and to 

maintain the ordinal nature of the data. The equations for volume growth of a structure in relation to Tanner stage are: 


Linearmodel:Volume=Intercept+f1TannerQuadraticmodel:Volume=Intercept+$1 Tanner+B2Tanner2Cubicmodel:Volume=Intercept+B1Tanner+B 
2Tanner2+P3Tanner3 


where 6}, Bz and B3 represent the constant terms defining the effects of each fixed term. To determine whether a cubic, quadratic or linear 
growth model with Tanner stage best fit the sample, an F test was performed on models where the marginal p-value of the highest order vari- 
able was significant (p < 0.05). 


Given that pubertal development is a developmental process, Tanner stage is necessarily highly correlated with chronological age, and this was 
found to be the case in our sample (r = 0.89 for males and r = 0.88 for females). It was therefore necessary to consider the effects of age on the 
developmental trajectories of the subcortical volumes, to see whether this could explain all of the puberty-related effects, and to establish 
whether models including age with puberty were a better fit than puberty models alone. Therefore, age was incorporated into the model to as- 
certain whether there were (i) dissociable main effects of Tanner stage and age, and/or (ii) a Tanner stage by age interactive effect. The model 
including the main effects of both Tanner stage and age and the tanner by age interaction is referred to as the combined model 

(Volume = Intercept + 8;(Tanner) + Bz (Age) + B3 (Tanner * Age)), whilst the model including only the Tanner stage by age interaction is re- 
ferred to as the interaction-only model (Volume = B; (Tanner * Age). Lastly, an age-only model was estimated (using linear, quadratic and cubic 
growth options as above). 


Comparison of Tanner-only, combined, interaction-only and age-only models was undertaken using likelihood ratio tests where possible. Where 
LR tests could not be performed, ie. if the models were not nested, Akaike Information Criterion (AIC) values were used to compare models. 
Models were considered to be a significantly improved fit if the difference between AIC values was 5.9 or greater, equating to an Akaike weight 
of the poorer model of less than 0.05. If the difference exceeded this, the inferior model was discarded, leaving only the models that were 


equally valid based on relative likelihood (Wagenmakers and Farrell, 2004). 


Results 


Demographic details for all participants are included in Table 1 and Fig, 1. The raw data for each sex and each structure are displayed in Fig. $2 
(Supporting Information). 


Models using Tanner stage as an explanatory variable for subcortical volume change 


Pubertal development, as measured by Tanner stage, was significantly related to the structural development of all six subcortical regions in both 
males and females (See Table 2 and Fig, 2). In both sexes, amygdala and hippocampus volume increased across puberty, whilst the other struc- 
tural volumes (NA, caudate, putamen and GP) decreased (See Fig. 2). In females, the growth trajectories were either linear (NA, GP) or quadratic 
(amygdala, hippocampus, caudate and putamen), whilst in males the trajectories were linear (hippocampus, GP), quadratic (putamen) or cubic 
(amygdala, NA, caudate) (See Table 2). The proportional volume change over puberty varied between structures from a 7.5% increase in male 
amygdala volume, to a 9.9% reduction in female GP volume, with some regions showing more modest volume changes e.g. caudate 2.3% reduc- 
tion in males. See Table 2 and Fig, 2 for full results. 


Models using Tanner stage and chronological age as explanatory variables for subcortical volume change 


For the amygdala, hippocampus and putamen in both sexes, as well as the caudate in females, a combined model (including Tanner stage and 
age as main effects, and a Tanner stage by age interaction) provided a significantly better fit than the Tanner-only model (See Table 3). For each 
of these structures, individuals in later stages of puberty (i.e. higher Tanner stage score) than their age-matched peers had volumes that are fur- 
ther along the developmental trajectory than peers in earlier stages of puberty. For females, individuals who were more pubertally mature than 
their age-matched peers had a larger volume in these structures in late childhood and early adolescence, an earlier peak volume and then a 
smaller volume until the end of puberty than their less mature, age-matched peers (see Fig. 3). In males, for the amygdala and hippocampus, in- 
dividuals in later stages of puberty (i.e. higher Tanner stage score than age-matched peers) had larger structural volumes throughout our stud- 
ied age range (7-20 years) than age-matched peers in earlier stages of puberty (see Fig. 3). In contrast, for the putamen in males, individuals in 
later stages of puberty had smaller structural volumes throughout our studied age range than age-matched peers in earlier stages of puberty 
(see Fig, 3). 


For the NA and GP in both sexes, and for the caudate in males, an interaction-only model provided a significant fit to the data, but the model was 
no improvement over the Tanner stage only model (See Table 3). Based on this interaction model, individuals in a later stage of puberty than 
their age-matched peers had a smaller NA or GP (both sexes) or caudate (males only) volume throughout the age range investigated (7— 

20 years) (Fig.3). 


Models using chronological age as an explanatory variable for subcortical volume change 


For some structures in males, the hippocampus and the GP, age-only models gave a significantly better fit of the data than models incorporating 
Tanner stage (hippocampus LR test compared to interactive model = 4.69, p = 0.030; GP LR test compared to interactive model = 9.29, 

p = 0.002). There was no statistically significant age-only model for volume change in the caudate in males over our age range. For the remain- 
ing structures in males (amygdala, NA and putamen) and all the six investigated structures in females, an age-only model did not improve the 
model fit over the models incorporating puberty measures. 


Discussion 


In the current study, we examined the relationship between puberty and growth trajectories of subcortical regions during adolescence using a 
large longitudinal dataset. Each of the regions studied showed significant longitudinal associations with puberty (Fig, 2). For many of the struc- 
tures (Females: amygdala, hippocampus, caudate, putamen; Males: amygdala, putamen), models including both age and Tanner stage modelled 
development better than Tanner-only models. This suggests that both variables are important factors when modelling volume change over ado- 
lescent development. Some structures in males (hippocampus, GP) were best modelled using only age as an explanatory variable, whilst the cau- 
date in males was best modelled using only Tanner stage. 


Despite the close proximity of the subcortical structures explored, there were clear differences in their structural development during adoles- 
cence. In both males and females, the amygdala and hippocampus continued to increase in volume during puberty, whilst the other structures 
examined decreased in volume, with structures changing between 2.3% and 9.9% across adolescence (Fig. 2). These results may reflect the dif- 
ferent mechanisms that influence macroscopic volume changes between structures, Alternatively, the regions may undergo similar growth pat- 
terns, but do so at different chronological time-points, resulting in different growth patterns within our restricted age and developmental range. 


Observation of pubertal development during adolescence provides an indirect marker of an individual's systemic sex steroid hormone expo- 
sure. The principal hormones involved are testosterone and dehydroepiandrosterone (DHEA), both androgens, and oestradiol, an oestrogen. 
Rising hormone levels, triggered by re-activation of the hypothalamic-pituitary-gonadal axis, leads to the series of physical changes classically 
associated with puberty. Androgens signal the development of adult-type body hair and skin changes in both females and males, in addition to 
gonadal development in males, whilst oestradiol primarily affects females, and causes breast and gonadal development. Tanner staging cate- 
gorises pubertal maturation by describing five stages of development from pre-pubertal to full maturation for each of pubic hair development, 
genital development and breast development. An individual's puberty stage is related both to how long they have been exposed to the sex 
steroid hormones, and to their current level of hormones (Dorn and Biro, 2011; Shirtcliff et al, 2009). Our results, showing that growth trajecto- 
ries of subcortical structures were related to pubertal development, suggest that these same pubertal hormones influence structural brain 
growth. 


Systemic pubertal hormones cross the blood-brain barrier in small concentrations (Marynick et al., 1976), and systemic concentrations of 
testosterone have been shown to be related to amygdala volume in both males and females (Neufang et al., 2009). Both androgens and oestro- 
gens induce synaptogenesis and synaptic pruning in rat and non-human primates (Ahmed et al., 2008; Hajszan et al., 2008; Sato et al, 2008) and 
it is likely that this process also occurs in humans, modulating brain growth across puberty (Matsumoto, 1991). Sex hormone receptors for both 
oestrogens and androgens are found throughout the brain in varying concentrations, with high levels in subcortical regions, particularly the hip- 
pocampus and amygdala (Abdelgadir et al., 1999; Clark et al, 1988; Shol] and Kim, 1989). The difference in regional receptor concentrations in 
the brain potentially helps to explain why different patterns of growth are seen across structures, and the resultant sexual dimorphism reported 
to emerge during adolescence in some regions (Brain Development Cooperative Group, 2012; Lenroot et al, 2007; Neufang et al., 2009; Sowell 
etal, 2002). 


We found that both males and females demonstrated pubertal effects related to amygdala growth, with increases in volume over puberty in 
both sexes. Although the overall volume change was similar between the sexes, the growth trajectories were quite different (Fig. 2), with females 
showing a large increase in volume in early puberty before peaking and decreasing, and males showing an increasing volume until the end of 
puberty. This is consistent with the different patterns of testosterone concentration in puberty, with males showing larger increases in concen- 
tration and a longer period of increase, and females showing a much smaller rise and earlier plateau of testosterone concentration (Ankarberg 
and Norjavaara, 1999). The smaller increase and earlier peak in testosterone concentration in females compared with males might explain the 
differences in trajectories between the two sexes, where the changes in neural structure that are seen are being modulated by systemic testos- 
terone concentrations (Zuloaga et al. 2008). This connection might reflect a direct effect of testosterone on amygdala volume, via the testos- 
terone receptors found in high concentrations in the amygdala, or indirect effects e.g. via aromatisation and effects on oestrogen receptors 
(Schwarz and McCarthy, 2008), or through interaction with growth hormone (Meinhardt and Ho, 2006) and its receptors present in the brain. 
Previous cross-sectional studies investigating changes in amygdala volume with puberty have reported volume increases with increasing puber- 
tal stage (and testosterone level) in males (Bramen et al,, 2011; Neufang et al, 2009), but have had conflicting results in females, where both in- 
creases (Neufang et al, 2009) and decreases (Bramen et al, 2011) in amygdala volume have been reported. The different findings between 
these two papers may result from the relatively small sample sizes or the differing age ranges of the two studies (Neufang et al., 2009: n = 46, 23 
male, age 8-15 years; Bramen et al, 2011: n = 80, 32 male, ages 10.7-14.0 years). As our sample shows, there is large variation in brain volumes 
between participants, and small samples may not be able to reflect the variation in the population, and within each group (early vs late puberty 
or Tanner stages 1-5). Our results indicate that both age and puberty impact on the structural development of the amygdala in males and fe- 
males, and differences in the findings of these two previous studies may reflect that they were sampling participants of different ages. 


In the hippocampus, the best model to explain growth in females was the combined model incorporating puberty and age, whilst in males the 
age-only model provided a better fit. This is consistent with an oestrogen-modulated growth pattern. During puberty, oestrogen concentrations 
in females increase by 4-9 times (Ikegami et al, 2001), and oestrogens are dominant in many pubertal changes. In contrast among males, rela- 
tively minor rises in oestrogens are seen during puberty due to aromatization of testosterone. The hippocampus in non-human primates has 
high levels of oestrogen receptors (Sho]] and Kim, 1989), and the increasing concentration of oestrogen that occurs predominantly in females 
coincides with the growth trajectories seen in the hippocampus in this study. 


In our analysis, both males and females showed a decrease in NA volume with puberty (see Fig. 2), and both sexes show on-going development 
throughout adolescence. Both androgens and oestrogens modulate the function of the NA, changing levels of dopamine release (Thompson and 
Moss,1994), and these results are consistent with the existence of macroscopic structural effects of both hormones on NA volume. For the NA, 
there was no statistical difference between the amount of volume change explained by the model based on puberty status alone, the model 
based on chronological age alone, or the model incorporating the interactive effects of both puberty and age. This may relate to the high corre- 
lation between age and pubertal stage. Note that such collinearity reduces the precision of the estimates of coefficients but does not affect esti- 
mates of model fit. Importantly, all the models appear to show on-going structural changes across adolescence in the NA. 


The caudate, putamen, and GP lie in close proximity to one another, and have related functions. Growth trajectories of these three regions 
across puberty are similar, with decreases in volume seen with greater pubertal development. The relative changes in volume for the caudate 
and putamen are the smallest of all the structures over the course of puberty (Fig. 2). The caudate in males was the only structure analysed that 
showed no significant relationship with age, emphasizing the importance of considering alternative variables that influence development over 
this age range. Previous cross-sectional studies of the caudate have shown similar findings, with weak or no age-related changes (Brain 
Development Cooperative Group,.2012; Ostby etal, 2009; Sowell et al. 2002); our results using a longitudinal dataset support this. Less is 
known about the sex hormone receptor concentrations in the caudate, putamen and GP. Sex hormones have been shown to impact upon recep- 
tor densities in these regions (Sumner and Fink, 1998), giving one potential indirect mechanism for the changes in structure seen, but further 
work exploring how sex hormones influence macroscopic structural changes is needed to help ascertain this. 


One of the major strengths of this study is the longitudinal nature of the dataset, which enabled us to model growth trajectories based on the 
trajectories of real individuals instead of using cross-sectional points. There is wide variability in structural volumes in the brain between indi- 
viduals, and repeated measurements of the same individual will therefore produce more accurate trajectories than assuming that cross-sec- 
tional data can be extrapolated to define trajectories, making the analysis more powerful to detect small but significant changes in brain volume. 
The structural development of the brain is likely to be affected by a number of variables and the use of a large longitudinal dataset maximises 
our ability to characterise the relationships of different variables with changing structural volume. The longitudinal nature of the data, the repli- 
cation of previous cross-sectional findings and the presence of plausible biological mechanisms allow us to hypothesize that puberty may have 
causal effects on brain development. However our study does not allow us to address whether the effects of puberty on structural brain growth 
are direct, via hormonal action on the brain tissue, or indirect, with the brain structure being shaped by how young people undergoing puberty 
may be treated differently in society than their less well-developed peers, or how they might perceive themselves differently (Blakemore etal, 
2010). To our knowledge, previous research has not tackled this subject, and further work exploring how pubertal hormones influence brain 
structural development in both animal and human studies, and combining hormonal and pubertal stage variables, may help to establish how 
this relationship is modulated. 


Our data are subject to a number of limitations. Self-reported Tanner stage, whilst the most widely used and clinically validated measure of pu- 
bertal development, is nevertheless a relatively crude measure which subjectively categorises puberty into broad developmental stages. It is 
therefore limited in its capacity to document accurately small developmental changes, and has significantly poorer resolution than age as an ex- 
planatory variable. This may mean that our puberty-related effects on growth trajectories are in fact underestimates, and further studies may 
help to validate our results and explore further their implications. Despite the large sample size, and longitudinal nature of the current dataset, 
there are limited numbers of participants at the extreme ages for each Tanner stage. This would be expected because our dataset is representa- 
tive of a normal population undergoing typical development, with a normal distribution of ages for each puberty stage, and therefore relatively 
small numbers of individuals at the extreme ends of normal puberty timing. These smaller numbers may reduce the accuracy of the model at 
these extreme ages and further research targeting narrower age ranges, and focusing on these extremes of normal pubertal development, and 
using different methods to measure pubertal maturation, would help to validate and further expand on our findings. Polynomial models, as used 
in this study, have been shown to be susceptible to the age range and age-centring used (Fjell et al. 2010) when performing age-based analyses 
of development. For this reason, we developed a clear a priori method for our analysis, based on our primary aim to explore the relationship 
between puberty and structural development of subcortical regions. Our age range was clearly defined, based on the largest available recent US- 
population based studies of pubertal timing (see page 6) to incorporate the reasonable ages associated with pubertal development. We included 
the full range of pubertal variation, our primary variable of interest, in the analysis, which should reduce any impact on the reliability of the 
model fit. A further limitation is our use of automated segmentation software to extract structural volumes for our six regions of interest. This 
method was chosen in view of the large number of scans (711) included in the study, and is widely accepted to be appropriate for very large 
scale-studies where manual tracing techniques are prohibitively time-consuming and resource intensive. Correlations between amygdala and 
hippocampus volumes for Freesurfer vs. manual tracing are high (Fisch] et al. 2002; Morey et al. 2009; Tae et al,, 2008). 


Conclusion 


We have shown that the structural development of subcortical brain regions is related to pubertal development during adolescence. This rela- 
tionship likely reflects the effects of systemic sex hormones on structural brain development. Examining brain development in relation to puber- 
tal development may provide additional information regarding the control mechanisms behind adolescent brain development, and in particular 
may shed light on how many of the behaviours classically associated with puberty come to arise. It may also help explain the development of a 
marked sexual dimorphism in psychiatric disorders around the time of puberty (Zahn-Waxler et al,, 2008), as there is emerging evidence of as- 
sociations between volumes of subcortical structures and psychiatric diagnoses (Karchemskiy et al., 2011; Rigucci et al, 2010). 
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Showing the age and Tanner stage of each participant at each study time-point. Males are shown in turquoise, and females are shown in orange. 
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Showing the growth trajectories for each subcortical region modelled against Tanner stage in females and males. For both sexes, the amygdala and hippocampus increase in volume 
over puberty, whilst the NA, caudate, putamen and GP decrease in volume. Shows the models for each region separately. Pink lines represent females and blue lines represent males. 
The solid line represents the best model fit, with 95% confidence intervals shown by the dashed lines. 


Fig. 2B. Showing the growth trajectories for each subcortical region modelled against Tanner stage in females and males. For both sexes, the amygdala and hippocampus increase in 
volume over puberty, whilst the NA, caudate, putamen and GP decrease in volume. Shows the models in terms of % volume change to allow comparison between structures. For each 
structure, the percentage volume was calculated for each pubertal stage as a proportion of prepubertal volume (at Tanner stage 1). 
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Showing how subcortical volumes change with age and puberty stage in (A) females and (B) males using the best-fit combined or interactive model for each structure. 


6000 


Age is presented on the x-axis, whilst puberty stage is indicated by the coloured lines (Red — Tanner stage 1, Yellow — Tanner stage 2, Green — Tanner stage 3, Blue — Tanner stage 
4, Purple — Tanner stage 5). Data for each graph were extracted from the combined models by calculating Intercept + (coefficient for main effect of Tanner x Tanner 

stage) + (coefficient for main effect of Age-centred x Age-centred) + (coefficient for interactive effect of Tanner-by-age x Tanner stage x Age). Age ranges for each Tanner stage line were 
decided using the ages and pubertal variation in our sample (see Fig. 1 for range). The age-only model for each structure is included (black dotted line) to aid interpretation. 


*For the caudate in males, there is no significant mode! that explains the developmental trajectory. 


Table 1 
Showing participants’ demographics. 


*This shows the number of participants who underwent an MRI scan at each Tanner stage, and the total number of scans collected of participants at each Tanner stage. Since some 
participants had more than one scan at a single Tanner stage, there are more scans than participants. 


Tanner stage _ 
Females Total sample 
N subjects/scans 


113 (12.0) 
41(172) 4198.2) 


Table 2 


Showing the Tanner-only best fit model for each of the six subcortical regions in females and males with the F test and p-value of the highest order variable. For each region, volume 
change across puberty is given in absolute (mm?) and relative (% change for the structure compared to initial volume) terms. 


Best-fitting model Volume change | across puberty Significance of highest order variable p-value 
Absolute (mm*) = % change 


Females 

Amygdala Quadratic 80 5.3% Fa,177) = 6.38 0.012 
Hippocampus Quadratic 191 5.1% Fu,177 = 6.26 0.013 
Nucleus accumbens Linear -50 -8.6% Fo, 178) = 18.73 < 0.0001 
Caudate Quadratic -131 -3.6% = Fas =10.78 0.001 
Putamen Quadratic ~ 282 -5.5% Fu, 7 = 16.15 0.0001 
Globus pallidus Linear -149 - 9.9% Fy, 178) = 45.02 < 0.0001 
Males 

Amygdala Cubic 130 75% Fea,2sa) = 7.40 0.007 
Hippocampus Linear 296 6.3% Fea, 256) = 69.01 < 0.0001 
Nucleus accumbens Cubic -55 -8.7% — Fra,259) = 4.44 0.036 
Caudate Cubic -92 - 23% Foy2s4) = 4.21 0.041 
Putamen Quadratic 159 -2.6% = Fu,255)= 7.91 0.005 


Globus pallidus Linear =79 -4.7% Fe, 256) = 19.22 < 0.0001 


Table 3 


Showing best fit models using a combination of Tanner stage and chronological age variables, with likelihood ratios and differences in AIC, Structures in BOLD show structures 
where the mixed Tanner stage and age model is a significantly better fit than the Tanner stage only model. 


a: Combined model refers to a model incorporating independent effects of Tanner stage and chronological age as well as an interactive Tanner stage by age effect. Interactive model 
refers to a model using only an interactive Tanner stage by age effect. 


b: In females, the best fit model for the putamen was a combined model including independent effects of Tanner stage, linear and quadratic chronological age as well as an interactive 
Tanner stage by age effect. 


* For these structures, a likelihood ratio test is not valid as the models are not nested and have the same number of degrees of freedom. Therefore, significance of the models has 
been judged using AIC differences. If AIC difference 2 5.9, the model is a significantly better fit (equivalent to Akaike weight of < 0.05). 


Structure Best-fitting Tanner and age model* Likelihood ratio test compared to Tanner only model p-value Difference between AIC 
Females 

Amygdala Combined 7.67 0.006 

Hippocampus Combined 23.96 < 0.0001 

Nucleus accumbens Interactive 3.21* i 3.21 
Caudate Combined 14.16 0.0002 

Putamen Combined” 28.13 < 0.0001 

Globus pallidus Interactive 3.43" “ 3.43 
Males 

Amygdala Combined 9.48" id 9.48 
Hippocampus Combined 23.06 < 0.0001 

Nucleus accumbens Interactive 4.63 0,099 

Caudate Interactive 5.29 0.071 

Putamen Combined 5.16 0.023 


Globus pallidus Interactive 1.51* * 1,51 


